Faraday rotation effect in periodic graphene structure 
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We report the magneto-optical rotation effect in a periodic graphene-sheet structure. Due 
to the masslessness of carriers in graphene, the magnetic response is very sensitive and the 
magneto-optical rotation effect is therefore significant. We predict that the Verdet constant 
of the periodic graphene-sheet structure is roughly 10 — 100 times that of rare-earth-doped 
magneto-optical glass in the infrared region. 



1 INTRODUCTION 

In the early 1960s, magneto-optical (MO) properties^ 
of materials started to develop rapidly because of the 
work of Bell Laboratories. Theoretical improvement 
and a lot of MO material synthesis led to developing 
many MO devices, such as MO modulators, MO iso- 
lators, MO sensors, magnetic optical circulators and 
MO memory, and even magnet ometryP However, ro- 
tatory power is not very high in usual materials. For 
instance, even for typical high- Verdet materials, such 
as MO glasses, the Verdet constant (VC) is only on the 
order of 0.1 — 1 min/ (Oe • cm)P Realizing a large VC 
is still a theoretical and experimental challenge. 

Surprisingly, Ref. 4 finds that we can obtain sev- 
eral degrees of rotation of the polarized direction when 
a linear electromagnetic wave passes through a sin- 
gle atomic layer of carbon, graphene. Recalling that 
graphene is the ultimately thin material in condensed 
matter physics, we regard this as a significant effect. 
There are also other studies^"^ ^ on this topic. 
For instance, Ref. [5] studies the phenomenon using the 
equation of motion, and Ref. [6] shows that graphene 
exhibits extremely broadband nonreciprocal polariza- 
tion rotation at subterahertz frequencies. 

So far studies have focused on rotation in mono- 
layer, bilayer or multilayer graphene. To make a prac- 
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tical MO device, that is, to obtain a rotation angle 
on the order of one radian even when the modulating 
magnetic field is not strong, one may resort to bulk 
graphite, for instance, flexible graphite. 



Here we study Faraday rotation from a different 
aspect. We first construct a periodic structure using 
graphene sheets. Assuming that the periodic length 
is far smaller than the wavelength of incident light in 
the structure, the periodic graphene structure (PGS) 
can be considered a medium. In this paper we study 
Faraday rotation when an incident linearly polarized 
electromagnetic wave passes through the medium (or 
PGS). 



In such a PGS a collective resonance mode, a plas- 
mon excitation, has a role. While there are only two 
parameters in monolayer or multilayer graphene,^ in- 
cident wave (angular) frequency uj and cyclotron fre- 
quency uo Cl there is a third parameter in bulk graphene, 
plasmon frequency uo pi which complicates our analysis. 
However, we will show that one can still produce an 
MO device with a high VC. 
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2 ANALYSIS OF FARADAY 
ROTATION EFFECT 

In semiclassical theory the VC is^ 
2m e c duo 

where c, /i, e, m e , n and uj are the velocity of light in 
vacuum, permeability of the medium, charge, electron 
mass, refractive index of the medium, and angular fre- 
quency of the incident light respectively. The equation 
tells us that to increase the VC, one may choose strong 
dispersion materials, for instance, MO glass. However, 
noting that ^ (multiplied by Bq) is the cyclotron an- 
gular frequency, we have another approach to increase 
the VC. 

As shown by Eq. (pQ), the strength of the MO rota- 
tion effect in materials is determined by the difference 
in velocity between the right circularly polarized and 
left circularly polarized waves. This difference reflects 
the asymmetry of the carrier response to an external 
magnetic field. In the medium, the carrier response to 
the external magnetic field is described by cyclotron 
motion, uj c = , where Bo and m e ff are the ex- 
ternal magnetic field and carrier effective mass respec- 
tively. However, from the expression for uo Cl one can 
also reduce carrier effective mass to enhance the car- 
rier response to the external magnetic field. 




Figure 1: Periodic structure of graphene 
sheets. The sides of the block cell are bi, 
b2 and d, which are parallel to the x, y and 
z axes, respectively. We impose the modulat- 
ing magnetic induction Bo and propagation 
direction of the electromagnetic wave both 
along z-axis. 



In usual materials m e ff ~ ra e , and such enhance- 
ment is not significant. Notice the massless Dirac 
fermion behavior of quasielectrons in graphene J^ i.e. 
m e ff = 0, which means the response to the external 
magnetic field should be sensitive. It is then natural 
to use graphene to get a high VC. In this section we 
study the Faraday MO effect in a PGS and find that 
we can use this PGS to get a very high VC. 

We first construct the PGS as in Fig. 1. The ele- 
mentary cell for such a PGS is defined as a rectangular 
block. The sides of the block are bi, b2 and d, which 
are parallel to the x, y and z axes, respectively, and all 
of them are far lower than the wavelength in the PGS. 
The graphene sheet outspreads in the x — y plane in 
each elementary cell. The ratio of the area of graphene 
to the area of a primitive cell in the x — y plane, 6162, 
is r. Here we just consider the case of r = 1, that is, 
each graphene sheet completely covers the x — y plane. 

There are two types of fields in the structure 
(medium). One is the modulating static magnetic field 
Bo, and the other is the electromagnetic wave field with 
E and B. Here we assume that the modulating mag- 
netic fields are far stronger than the electromagnetic 
wave fields, i.e. Bo ^> B. We impose the modulat- 
ing magnetic induction Bo and propagation direction 
of electromagnetic wave both along the z-axis. Hence 
B and electric field E are in the x — y plane. 

Landau diamagnetism is negligible because there is 
no external electric field. Therefore, we approximate 
the relative permeability as \i r ~ 1. 

Now we use the concept of an effective plasmon to 
study dielectric permittivity. We first write the effec- 
tive carrier density, n = n2r/d, where r\2 = £^ (^7) 2 
is the number of carriers per unit area, Ep is the Fermi 
energy (or chemical potential at zero temperature) 
without external fields, vp = 1 x 10 6 m/s is the Fermi 
velocity of graphene and g = 4 is the degeneracy!^ In 
the graphene sheet, the Fermi energy Ep and Fermi 
momentum pp satisfy pp = Ep/vp = h^/^n^- 

Notice that since the Fermi velocity vp = <C c, 
the long-wave limit is always correct in the periodic 
structure. If we turn off the modulating field, that is, 
we ignore the cyclotron motion, the dielectric permit- 
tivity can be written as a scalar JED EH ESI 

e(«) = l-^, (2) 

in the long- wave limit, where u p is the plasmon angular 
frequency or resonant angular frequency. We empha- 
size that, as shown in Refs. [131 [14] and [T7J this form 
for dielectric permittivity is valid both in the quan- 
tum and classical theories, but the expression for u p in 
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the quantum theory is different from that in the clas- 
sical. In the classical theory, u p - 

graphenef^ 1 ^^ u p = yj^^f- due to quantum and 
'relativistic' effects, where eo is the vacuum permittiv- 
ity. In this case, the permittivity is isotropic the in 
x — y plane and there is no rotation effect. 

Now turn on Bo along the z— axis and discuss the 
rotation effect. To make semiclassical theory apply, 
we first assume uj c = ^^,uj p ,uj <C ^f. 16 Since 
fthj <C Ep, the photon energy is much lower than 
Ep, so we consider only the cyclotron motion near the 
Fermi energy. This is because the carriers with en- 
ergy far lower than the Fermi surface cannot absorb 
energy from a photon because of the Pauli exclusion 
principle. Under the assumption, as pointed out by 
the Appleton-Hartree theory, we write the relative di- 
electric permittivity e r in the x — y plane as^ 



(3) 



where e\ 



= 1 

eB 



^2 
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uj(uj2-uj 2 ) 



eB v F 
Pf 



We have shown in Ref. fTTl 



and uo c 

that u p and u c are quantum angular frequencies of 
the plasmon and cyclotron respectively. Equation (|3j) 
tells us that the effective permittivity is determined by 
two physical behaviors in the PGS, namely, cyclotron 
motion and plasmon excitation, characterized by two 
frequencies, u c and uo p respectively. The cyclotron mo- 
tion is determined by the modulating magnetic field 
and the plasmon is a collective mode. On one hand 
we find that, if uj c <C co>,co> p , the dielectric permittivity 
tensor will degenerate into Eq. (|2]), while on the other 
hand, e 2 ex uo c /uo will lead to the optical rotation effect. 
Maxwell's equations in the PGS are 
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and constitutive relations are 

/i H, 
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where /io is the vacuum magnetic permeability. 

Supposing the electric field and magnetic field can 
be written as B(E) = B (E )e iu;£ - i ^, where B, E, B 
and Eo are vectors in the x — y plane, one finds that 
for the electric field 
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where the identity fio^o = 1/c 2 is used. The relation 
for the magnetic field is similar. For each u, cor- 
responding to solutions of nonzero E = (E x , E y ) in 
the above equation, there are two eigenmodes: k\ = 
^V e i + e 2 with eigenvector Eoi = (e x — \e y )A/2 and 
&2 = ^ V 6 i — e 2 with eigenvector E02 = (e x +ie y )A/2, 
where A is a constant. The former eigenmode is left 
circularly polarized while the latter eigenmode is right 
circularly polarized. Since k\ 7^ hi for a certain gj, 
the velocities of these two circularly polarized waves 
are different. Therefore, a rotation effect occurs for a 
linear electromagnetic wave passing through the PGS; 
that is, if the incident wave is linearly polarized from 
the negative z-axis, the emitting wave is a still linearly 
polarized but with a rotational polarized direction. 

To analyze the rotation effect, we first suppose that 
both left and right circularly polarized waves can prop- 
agate in the medium and write out the electric fields 
of these waves as 



Ei |-(cos(cjt 
E 2 = |(cos(cjt 



k\z)e x + sin(cjt — k\z)e y \ 
k2z)e x — sm(ujt — k2z)e y \ 



(7) 



with total E = Ei + E 2 . This expression means that 
in the z = plane the linear polarization is along the 
x— axis. We get 



E 



Acos(ujt — 
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Therefore, the total wave vector is 



k = — (Vei + e 2 
Zc 

and the rotatory power is 



e 2 J 



£2 - Vei + £2)- 



(8) 



(9) 



(10) 



To consider rotatory power in detail we introduce 
two dimensionless quantities, Q = uo c /uo p and ujo = 

We 



_Q_ 



o; (Q 2 -^o) ' 



uj/uj p . Then, e\ = 1 + n2 _ 2 and e 2 

find that when uj > ^(Q + ^Q 2 + 4), both e 1 ±e 2 > 0, 
and a rotation effect occurs. Notice that if we take 
k < 0, the rotatory power is still a. This phenomenon 
is the well-known non-reciprocity of the MO effect. 
Now, we have 



ca 





Q 
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In fig. 2(a) we show vs. uo with different Q. We 
find that for fixed (high) uo and lj p , \a\ increases with 
decreasing ujo or increasing Q (and thus Bo). The 
behavior agrees with that of multilayer graphene, as 
Fig. 3 in Ref. 4 shows. (The relative negative sign 
is due to the fact that we take ep > while in Ref. 

uJ 2 — 1 

[U £f < 0). At Q = , a reaches its maximum 

value, ^ = ujq^J -^jt^, which can also be seen in Fig. 
3. However, unlike single- and multilayer graphene, 
one cannot get a further step at this time. If we fur- 
ther reduce the incident frequency, one of the circularly 
polarized waves should decay in the medium, and the 
emitting wave should generally be elliptically polarized 
if the thickness of the medium is finite. Notice that for 
single and multilayer graphene, the decay is not impor- 
tant and we can also get a linearly polarized emitting 
waveP 




^0 COq 



Figure 2: ^ vs. ujq with different Q. 
(a) small Q and higher ujo. (b) larger Q 
and lower ujq. In figure (a) the dashed 
line, dotted line and solid line are for Q = 
0.01, 0.1, and 0.5 respectively and the red 

line is ujq J 4a; 2_ 2 . In figure (b) the dashed 
line, dotted line and solid line are for Q = 
1.5, 2.5, and 3.5 respectively and the red line 

Ordinarily, Q is a small quantity. Theoretically, one 
may also study the case of large Q. One finds that if 
Q > l/y/2 and ^( v / 4 + Q 2 - Q) < ujo < Q, both right 
and left circularly polarized wave can pass the PGS 
without decay. We show the result of lower frequency 
with larger Q at Fig. 2(b). Notice in this region rota- 
tory power a > 0. However, in this region the lower 
uj p and uj c make the quasielectron behavior dominated 
by quantization, thereby making our results doubtful. 

Figure 2(a) shows that \a\ is on the order of u p /c 
when the frequency of the incident wave is close to the 
resonant frequency. For instance, when we adjust the 
parameters to Ep = 0.3 eV, d = 205 nm and Q = 0.01 
(B ~ 4.7T), we find that Huj p = 0.09eV and \a\ is on 
the order of 10 — 10 2 r ad/mm, which is significant. 

Fig. 3 shows that at small Q, rotatory power is 
proportional to Q, or in other words, a oc Bq. From 



the Taylor expansion we get the VC with 

= |oj = ev F fi 

Bo/^o 2cE F uoy / ujQ - 1 
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= "T^TT^ - ^=(min/Oe- cm), (12) 

where we assume Ep = x(eV). In the structure we 
generally choose x ~ O.leV. When Q « 1, one way 
to improve V is setting The parameter to make fre- 
quency close to the resonant frequency, ljq ~ 1. If 
we set ujq^/uJq — 1 ~ 1 and set Ep — 0.3eV, V is on 
the order of 19min/Oe • cm. The VC of the graphene 
structure is about one or two orders larger than that 
of MO glass. This is thus a significant advancement in 
magnetic optics. 








Figure 3: ^ vs. Q with different ujq. The 
dashed line, dotted line and solid line are 
for ujo = 1.01, 1.05, andl.2 respectively. The 

red line is \{Q + ^ITO 2 "),/^ 3 . 

2V V Y 3Q+V4+Q2 

Equation (pQ) shows that V in a usual material is 
inversely proportional to m e ff, where m e ff is the ef- 
fective carrier mass. In graphene, m e ff = and the 
usual expression is not valid. However, Pf/vf plays 
the same role as m e f f , and thus one can obtain a very 
large (not infinite) VC. Here we estimate the order 
of enhancement. In usual materials, the carrier mass 
m e ff ~ m e c± 10 5 eV/c 2 . However, in graphene, cy- 
clotron motion mass plays the role of effective mass, 
m c = Ef/v f . If we choose Ep = 0.2eV, we find that 
m c - 10 4 eV/c 2 - 10- 1 - 10- 2 m e . Thus, the VC in 
the PGS is about one or two orders larger than in rare- 
earth-doped MO glass. 

Furthermore, as shown in Fig. 3, when Q tends to 

oj 2 — 1 

, the rotatory power exhibits nonlinear behavior 
and VC increases further. 

3 DISCUSSIONS 

In this manuscript we have discussed the MO rotation 
effect in a PGS. The masslessness of graphene carri- 
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ers gives them a very sensitive response to an external 
field, and the VC of the PGS is very high. We predict 
that the VC in the structure is about one or two orders 
larger than in rare-earth-doped magneto-optical glass. 

Furthermore, even if uj is not close to the reso- 
nant angular frequency uj Pi we may also improve B 
to increase the VC, although such improvement is con- 

uj 2 — 1 

strained by the condition Q < -° — . 

In a usual material, the effective mass is always 
a constant, whereas in graphene, the role of m e //, 
m c = Ep/v Fl is adjustable. This behavior can further 
expand the application range of the magneto-optical 
effect. 

We conclude that besides MO glass, we can also use 
a graphene-sheet periodic structure to attain a high 
VC, increasing the VC up to one or two orders of mag- 
nitude. We hope our study is helpful in designing new 
types of MO devices. 
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